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Abstract 



In this paper, we study the competition of linear and nonlinear lattices and its effects on the sta- 
bility and dynamics of bright solitary waves. We consider both lattices in a perturbative framework, 
whereby the technique of Hamiltonian perturbation theory can be used to obtain information about 
the existence of solutions, and the same approach, as well as eigenvalue count considerations, can be 
used to obtained detailed conditions about their linear stability. We find that the analytical results 
are in very good agreement with our numerical findings and can also be used to predict features of 
the dynamical evolution of such solutions. 

o . 

^ ■ 1 Introduction 



The field of Bose-Einstein condensates (BECs) in atomic physics [T] has provided a renewal of interest 
in solitary wave and nonlinear excitations. This is especially so because the leading order, mean-field 
description of interatomic interactions corresponds to a cubic nonlinearity. This, in turn, is responsible 
for the emergence and experimental observation of a variety of matter-wave solitons in BECs, including 
bright dark [I] and gap [5] solitons. Furthermore, atom optical devices have been proposed, such 

as, for instance, the atom chip [6] that would allow the possibility to controllably manipulate such 
nonlinear structures. 

It is well-known that bright (respectively, dark) matter-wave solitons arise in BECs with attractive 
(respectively, repulsive) interatomic interactions, i.e., for atomic species with negative (respectively, 
positive) scattering length a. 

One of the particularly appealing features of the BEC setting is the existence of a wide variety of 
experimental "knobs" that can be used to manipulate or control the relevant structures. In particular, 
interfering laser beams can be used to produce a standing wave pattern, known as the optical lattice, 
providing a periodic linear potential for the condensate. This type of structure offers a large variety of 
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interesting phenomenology including Bloch oscillations, Landau-Zener tunneling, dynamical instabili- 
ties, gap excitations among many others and a considerable amount of review works have already been 
dedicated to this topic O E] . 

On the other hand, magnetically-induced Feshbach resonances can be used to modify at will both the 
magnitude and the sign of the scattering length by tuning the external magnetic field; see e.g. |10| and 
also [3 13] where the Feshbach resonance in 'Li BECs was used for the formation of bright matter- wave 
solitons. The ability to modulate the scattering length has led to a large variety of studies where this 
mechanism has been used. For instance, time-dependent modulations of the scattering length were 
proposed as a means of preventing collapse in higher-dimensional BECs [11], or as a way of producing 
robust matter- wave breathers [12] . among others. 

A more recent suggestion has been to add to a constant bias magnetic field a gradient in the vicinity 
of a Feshbach resonance, allowing for a spatial variation of the scattering length, thereby providing what 
has come to be termed a "collisionally inhomogeneous environment". Notice that, given the availability 
of magnetic and optical (laser-) fields, the external trapping potential and the spatial variation of the 
scattering length can be adjusted independently (see [13] for more details on the relevant configuration). 
In this latter setting, a variety of propositions of interesting dynamical phenomena have been made 
concerning scenarios for the emission of solitons [H] , delocalizing transition of matter waves [15] , or the 
dynamics of the waves in random [IB] , linear [IT], periodic [18j[l9] or localized [20J spatial modulations. 
A number of more mathematically minded results on the existence and stability of waves have also 
appeared in [21] and a technique for analytically constructing exact solutions in [22J. 

The framework of collisionally inhomogeneous environments in combination with external optical 
lattices provides an ideal environment for competition. This was illustrated in examples of modula- 
tional instability of Bloch states [19] and of the delocalizing transition in one-dimension |15j . In the 
present work we develop these ideas exploring more general lattice profiles. In particular, external 
(linear) potentials and collisional (nonlinear) potentials in (as well as out of) phase will be consid- 
ered. It will be demonstrated that when in phase, these potentials provide a competition leading to 
a number of interesting effects including stabilization/destabilization thresholds and even the mutual 
annihilation of the two potentials to provide an effectively quasi-translationally-invariant environment. 
The effective potential landscape where the solitary waves (of the bright type) live will be obtained 
following the Lyapunov-Schmidt considerations of [23]. Then, the relevant (translational) eigenvalue of 
the linearization will be computed based on the curvature of this effective potential landscape and the 
stability /instability of the waves will be assessed (and the relevant transition points will be obtained). 
These results will be confirmed by a second independent method based on a direct count of unstable 
eigenvalues. Finally, we will examine the variation of a number of relevant key parameters (such as the 
amplitude, the wavenumbers or the relative phase) of the potentials, in order to evaluate the validity of 
the approach. Furthermore, when the structures are unstable, we will examine what this approach can 
suggest regarding the actual instability evolution dynamics. The understanding that we will develop 
will enable us to manipulate the ensuing solitary waves in such a complex territory and to understand 
their dynamical behavior in the presence of linear and nonlinear lattices. 

Our presentation will be structured as follows. In section II, we will present our analytical results; in 
section III, we will corroborate these results by means of numerical computations. Finally, in section 
IV, we will summarize our findings and present some interesting directions for future study. 
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2 Analytical Results: Solitary Wave Statics and Dynamics 



The prototypical framework in which we will consider the above discussed competition of linear and 
nonlinear lattices is that of the perturbed nonlinear Schrodinger equation of the form 

iu t = -^u xx - (1 + eni(x)) \u\ 2 u + en 2 (x)u. (2.1) 

In Eq. (|2.ip . (x,t) £ R x M + and u G C. While we will keep the presentation of the mathematical 
results as general as possible, the particular case of interest in the selection of the nonlinear and linear 
lattice will, respectively, be: 

ni(x) = Acosihx) . . 

ri2(x) = B cos(k2X + Acj)) . 

where A, B, k±, k 2 and Acj> are real constants. Notice that the lattices have the same functional form, 
which will allow us to reveal more lucidly the relevant competition between the corresponding terms. 

When e = 0, Eq. (|2.1I) has the well-known stable localized soliton solution given by 

u(x, t) = ^j2sedi[^jl(x - £)] e *M*-0+*] (2.3) 
where [i > 0, £ = vt is the position of the soliton center, v is the velocity of the soliton, and 5 = 

We presently focus on the stationary modes with v = 0. Given the monoparametric nature of 
the family of the respective soltuions, we can fix /x in what follows (in fact, we will fix [i = 2 in 
our numerical computations below). Because of the rotational and translational invariance of the 
unperturbed equation, this solution is unique only up to rotational and translational symmetry. 

On the other hand, when e > 0, the translational invariance of the equation is broken, which may 
naturally lead to the potential destabilization of the localized states, depending on the perturbation 
parameters. This is the problem that we will examine in what follows under the influence of both linear 
and nonlinear lattices. 



2.1 Hamiltonian Perturbation Approach 

The existence and nature of localized solutions to perturbed Hamiltonian systems, of which Eq. (12. lj) 
is a particular case, was studied in [23] (and subsequently in a broader setting in [24]). A general 
perturbative approach was developed in these works based on Lyapunov-Schmidt solvability conditions 
[25], and relevant stability calculations were formulated on the basis of the Evans function [26\ 127] . 
Here, we present some of the general features of the theory, adapt our problem to the general framework 
of [23j El] and subsequently apply these methods to the problem of interest. In order to apply these 
criteria, it is convenient to recast Eq. (|2.ip as 

du 5E , 

where E(u) = Eq{u) + eE\{u). Here, 

/+oo ^ 
- {\u x \ 2 - M 4 ) dx, (2.5) 
-oo ^ 
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and 



Ei(u) = j {n 2 {x)\u\ 2 - ^ni(x)M 4 ^ dx. (2.6) 

Then, for fixed [/,, the intuitive condition for the persistence of the wave is given by [23] 

V e £i(u) = 0, (2.7) 

where £ is the previously-free parameter associated with the invariance (in the case of translation, it is 
the center of the pulse (|2.3p ). This condition implies that the wave is going to persist only if centered 
at the parameter-selected extrema of the energy (which are now going to form, at best, a countably 
infinite set of solutions, as opposed to the one-parameter infinity of solutions previously allowed by the 
translational invariance). 

Equally importantly, from this expression and from the nature of the wave, one can infer stability 
information about the solution of interest. In particular, the stability of the perturbed wave is deter- 
mined by the location of the eigenvalues associated with the translational invariance; previously, the 
relevant eigenvalue pair was located at the origin A = of the spectral plane of eigenvalues A = A r + zAj. 
On the other hand, we expect the eigenvalues associated with the U(l) invariance (i.e., the phase in- 
variance associated with the L 2 conservation) to remain at the origin, given the preservation of the 
latter symmetry under the perturbations considered herein. To compute the relevant eigenvalues we 
refer to the framework put forth by the works of |23|. [24J (adapting the notation of the latter work) in 
the following form. Using Proposition 6.1 of [23], we expect that the perturbed system eigenvalues will 
be given by the matrix equation: 

det[Mi + \\D G } = (2.8) 

where 

(d x u°, -xu°) \ / ^ 1/2 

also, Mi is given by 



Dg= \ o ) = [ o V /2 ); (2,9) 



Mi = ( ' " I - | ' [I — ' ~ "FT! 3 "' " i }"■' " 1 (o 10) 

I o oil o 1 



In the formulation of [24] . the relevant eigenvalues are obtained to leading order as A = \/eAi, and u° 
denotes the solitary wave of Eq. (|2.3p . Hence, we conclude from Eq. (|2.8p that as indicated above, the 
eigenvalues associated with the rotational invariance will be preserved at A = 0, while the translational 
eigenvalue will be shifted according to: 

A 2 = ^pr [ ( - d ^(uy - l 7 ^(ur) dx (2.11) 

/xV2 J \2 dx 2 v ; 4 dx z y > J y ' 

Hence, the corresponding eigenvalue can be straightforwardly evaluated, provided that we first com- 
pute the extrema of the effective energy landscape E\, which, as a function of £, will hereafter be 
denoted as V e f /(£), or more precisely V e ff{€) = zE\- This V e ff will be the effective energy landscape 
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and the stability or instability of the configuration will be associated with the convexity or concavity 
of this effective energy landscape. 

We now proceed to evaluate the relevant expressions of the general theory for the special case of 
interest herein, namely for the potentials of Eq. (j2.2j) . The effective energy landscape can be evaluated 
after performing two straightforward contour integrations, that yield 

V 9f m = -exA* (kl + 4 M ) COS( ^ . + e,Bk 2 C ° S ^ + A ^ ( 2 .12) 

In the simple case of A = A(p = 0, the above result reduces to that in [23], which leads to the well-known 
conclusion of [28], according to which a maximum of a linear periodic potential leads to an unstable 
solitary wave configuration, while the opposite is true for a minimum of a periodic potential. However, 
in our case, there is an intriguing interplay between the -B-dependent term stemming from the linear 
optical lattice and the ^4-dependent term, emerging from the nonlinear optical lattice. This competition 
leads to the potential for stability-instability transitions for the wave, based on the properties of the 
trapping (such as (A, B, ki, k 2 , A<^>)), but also the properties of the wave itself (since the expression of 
(|2.12|) is explicitly dependent on ji). 



2.2 Eigenvalue Count Approach 

It is worthwhile to note, however, that the above stability results, based on the formulation of |23[ I2H] 
can also be alternatively derived using the approach of [28, 29J. We present this alternative formulation 
here, since we consider it to be a nice complement to the direct eigenvalue computation of [231 124"] using 
the Hamiltonian perturbation technique. It is well-known that the stability of the solitary wave [24j is 
determined by the number of negative eigenvalues of the operators 

(1 + eni)(u e ) 2 + en 2 + | (2.13) 

3(l+eni)(u £ ) 2 + en 2 + | (2.14) 

In particular, if n(L) denotes the count of negative eigenvalues and \n(L e + ) — n(L e _)\ > 1 [29], then 
the solitary wave is unstable. The superscript e in the operators is to distinguish the e ^ and the 
e = cases. In the latter, L+ is well known to have a single negative eigenvalue (— 3^/2) with an 
eigenfunction spanned by (u ) 2 and a single zero eigenvalue with an eigenfunction spanned by dvP/dx, 
while LP. has no negative eigenvalues and a single zero eigenvalue with an eigenvector spanned by u°. 
In the perturbed case, U_ retains its zero eigenvalue with eigenvector spanned by u e . Also, from the 
perturbation theory of Schrodinger operators [30| it is known that L e + has a negative eigenvalue near 
—3/j,/2, and a second eigenvalue near 0, which will be denoted by \l* . Both of these eigenvalues are 
analytic in e, at least in a neighborhood of the real axis. Then the stability question (given also that 
for this branch of solutions dN/dfJ, > |24j ) is completely settled by the count of negative eigenvalues 
of L e + . In particular: 



Li 



l_cP_ 
2dx^ 
ldP_ 
2dx^ 



If U + has two negative eigenvalues, then the solitary wave will be unstable since \n(L e + )—n(L e _ 
2; 

If LI has only one negative eigenvalue, then the coherent structure will be stable. 
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Hence, the stability issue hinges on the shift of the zero eigenvalue (corresponding to translational 
invariance, when e = 0) in the presence of the perturbation. Thus, similarly to [28J, we will consider 
the quantity (L e + u,u) with u being the eigenvector corresponding to the eigenvalue A^ ; here (v,w) 
denotes the L 2 -inner product of v with w. u can be decomposed asu = u\\ +u±, where ity is proportional 
to du £ /dx and (u±,du e /dx) = 0. Along the lines of (3TJ [32] one can show that \\u e — u°\\ H 2 — > 0, and 
therefore uu — > dvP/dx (up to a proportionality factor) as e — > 0. From perturbation theory [33J it is 
also known that the eigenvector u of L e + is analytic in e. Therefore u± — > 0, and actually u± ~ e, i.e., 
it will be of the order of the perturbation. Then, we have 

(L%u,u) = (L e +U]l , Ull ) + 2(L e +Ull ,u ± ) + (L e + u ± ,u ± ). (2.15) 

However, each of the second and third terms will be of order higher than the first (at least 0(e 2 ), while 
the dominant one will be of 0(e)), hence the solitary wave stability will be determined by (L^u»,u»). 
But then, 

(L>,u)=r 2 (L^,^) + 0( e 2 ); (2.16) 

r is an appropriate proportionality factor (between u\\ and du e /dx). By means of a direct computation 
(differentiating the equation satisfied by the stationary state) one has that 

TC du e dno c dn\ , . . 

which, in turn, forming the inner product with du € /dx, and integrating by parts leads to the key result, 
namely: 

(L%u,u)=rh J Q^(-°) 2 -^(-°) 4 )^ + 0( e 2 ) (2-18) 

The relevant integral term of the right hand side can be seen by direct inspection to be equivalent 
(up to a negative-definite proportionality factor) to the expression for A 2 in Eq. (|2.11j) . Its positivity 
(indicating a shift of the zero eigenvalue to positive values) will imply stability, while its negativity 
(indicating a shift of the zero eigenvalue to negative values) will lead to instability. This conclusion 
is fully equivalent to the ones obtained from Eq. (|2.1ip [although the latter, in some sense, contains 
additional information yielding a quantitative measure of the relevant eigenvalue] . 



2.3 Solitary Wave Dynamics 

In order to describe the dynamics of a soliton of Eq. (|2,ip at e <C 1 one can employ the perturbation 
theory for the NLS soliton [34j . or more precisely, the adiabatic approximation. If e ^ then £ ^ vt 
and has to be found from the equations of the adiabatic approximations. The straightforward algebra 
yields 

dH_ i dv eff (Q 

dt> ~ N di { ' yj 

where N = J \u\ 2 dx = is the number of particles [it is an integral of motion of (|2.1)) ]. 
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It follows from (|2.19|) that there exist different types of motion of the soliton. In particular if 
&2 = k\ = k and Acfi = the soliton dynamics reproduces the mathematical pendulum. The respective 
motion of a soliton can be either periodic or translational (i.e., unbounded), depending on the initial 
conditions. Another special case arises for A(k 2 + 4/i) = 12B, in which case the right hand side of 
(|2.19p becomes zero and in the adiabatic approximation £ = vt, i.e. the motion becomes linear because 
the periodic nonlinearity effectively compensates exactly the periodic potential. 

More sophisticated evolution scenarios can be observed for k2 ^ k\ depending whether they are 
commensurable or not. 

3 Numerical Results 

We now proceed to describe our numerical results comparing with the analytical prediction of the pre- 
vious section. We use e = 0.1 (for which we expect the perturbative description to still be meaningful), 
and vary the relative parameters of the two lattices (linear and nonlinear). 

Our first set of numerical results consists of setting B = 1 and k\ = k<i = 2ir/5, Ac/) = and varying 
A. In this way, we can test the validity of our predictions for amplitude variations (in this case of the 
nonlinear lattice). Our results are summarized in Fig. [TJ We have varied A £ [0, 2], finding that there 
is a stability change within this interval. In particular, the left panel of the figure shows the case of 
A = 0.5 which is unstable, and of the stable A = 2. The right panel shows this transition in terms of 
the real part (and also of the square) of the relevant eigenvalue associated with the translational mode. 
It is found that this eigenvalue pair starts out as real, for small A, and becomes imaginary for A > 1.33. 
We use Eq. (12. lip to theoretically predict this transition as occurring at A = 1.253. While we see 
that both qualitatively and fairly quantitatively the dependence of the eigenvalue on the parameter 
is captured accurately by our theoretical result, it is meaningful to rationalize the ~ 6% error in the 
critical point estimation. It is, in fact, observed that the soliton does not maintain its amplitude in this 
continuation process (as a function of ^4), but rather that its amplitude is reduced from ~ 1.455 for 
A = to 1.353 for A = 2. This clearly shifts the critical point upwards, whose analytical expression in 
this setting of k\ = k2 and Ac/) = can be easily seen to be A cr = \2B /{k 2 + 4^). This is in agreement 
(in fact, even quantitatively, if one uses the above amplitude variation) with what is observed in our 
numerical results. 

In the bottom panel of the figure, we show the result of the unstable dynamical evolution for the case 
of A = 0.5. It can be seen that as a result of the dynamical instability the solitary wave starts moving 
to the left, eventually executing oscillations between the two maxima of the effective potential of Eq. 
(|2.12p . In the same plot, we show the result of the adiabatic soliton perturbation theory in this case 
(this is a rather "stringent" test of the theory given the unstable dynamical evolution). We observe 
that the Eq. (I2.19P performs well in approximating the soliton trajectory over the first oscillatory 
cycle. However, for longer times, we observe it to gradually increasingly fail to capture the relevant 
oscillation. This can be seen to be due to the fact that the solitonic trajectory emits small wakes of 
radiation as it arrives at the turning points, resulting in a weakly damped oscillation, a feature which 
is not captured by our present considerations. However, we note in passing that methods similar to 
those developed by Soffer and Weinstein [35] can be used to rigorously account for such corrections. 

Our second parameter variation involved the role of the wavenumbers. In particular, for the results 
reported in Fig. [2J we have used A = B = 1, and fixed k<z = 2ir/5 and Acf) = 0, varying k%. One can see 
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Figure 1: The top left set of panels shows the cases of A = 0.5 (unstable, left) and A = 2 (stable, 
right). Both the solution profile (solid line) and the effective potential (dashed line) is shown in the top 
subplots, while the spectral plane of (A r , Aj) is shown in the bottom. The presence of a real eigenvalue 
indicates instability in the latter. The right panels show the "trajectory" of the relevant eigenvalue 
(real part in the top and squared eigenvalue in the bottom) as a function of A. The solid line indicates 
the numerical result while the dashed line the analytical prediction for the eigenvalue. The bottom 
panel shows the evolution of the unstable configuration for A = 0.5, in the effective potential (shown 
out of scale in the graph by a dash-dotted line). In this spatio-temporal contour plot of the squared 
modulus (in fact, its opposite is shown for clarity) of the solution, the result of the solitary dynamics 
of Eq. (|2.19p is superposed as a thick dashed line. It is clear that the solution, as a result of the 
instability, oscillates between two maxima of the effective potential. 
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Figure 2: Same as Fig. [H but for the case where k\ is varied (i.e., one of the wavenumbers, instead of 
the amplitude A). The top left panels show the solution (solid line) and effective potential (dashed line) 
for k\ = 1 (left subplots) and k\ = 2 (right subplots), as well as the corresponding spectral plane for the 
stability (bottom subplots). The top right panels show the theoretically predicted (dashed line) versus 
numerically obtained (solid line) eigenvalue of interest. The bottom panel shows the spatio-temporal 
evolution of the unstable case with k\ = 1, exhibiting oscillations in the effective potential landscape 
(shown out of scale by a dash-dotted line). 



that in this case, the effective potential landscape changes significantly in term of its local structure 
(in the previous example, it did not change, in that it was simply two cosinusoidal terms with different 
signs, so it was simply a matter of which had the largest "strength"). The other important feature is 
that in this case, as well, there is a transition from instability to stability, as k\ is increased. In fact, the 
theoretical prediction for the critical point is k\ = 1.433, while the numerical one is k\ = 1.46. Once 
again this can be seen for the two different settings of the left panel (the unstable case of k\ = 1 and the 
stable case of k\ = 2), and is captured extremely accurately by the prediction of Eq. (12.111) about the 
location of the relevant eigenvalue (associated with translation). Furthermore, once again, the effective 
potential landscape that can be computed from Eq. (|2.12p can provide very useful information not only 
about the stability of local extrema but also about the instability dynamics. The latter is observed in 
the bottom panel of the figure. Note, however, that while the effective potential predicts accurately 
the turning points of the solitary wave dynamics, the situation is more complicated with the dynamical 
equation of motion of (|2. 19j) . While, once again, the latter predicts very accurately the first oscillation 
cycle, its non-accounting of the radiative corrections of the motion leads to dynamics that overcomes 
the shallow potential barrier at x = 0; this is not true, however, for the full PDE dynamics. This 
should serve as a note of caution in regard to using the adiabatic approximation in such (marginal) 
cases. 
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Figure 3: Same as the top panels of the previous figures. The left panels show the configuration and 
its stability for Acj) = 2 (leftmost panel, the wave is shifted to the right) and A(p = 4 (where the wave 
is shifted to the left). The right panel shows the relevant unstable translational eigenvalue in the top 
subplot and the center location of the wave in the bottom subplot. The theoretical results are given 
by dashed lines and compare very favorably with the solid lines of the full numerical results. 



We also explored the role of the phase difference between the linear and nonlinear lattice, by varying 
A<p E [0, 2n], for A = B = 1 and k\ = k% = 2ir/5. One of the particularly interesting features of this 
example is that while the instability of the original configuration is not modified by this variation, the 
location of the solitary wave is. This is naturally expected on the basis of Eq. (|2.12p . In particular, 
we observe that the bright soliton's center location features an oscillation around £ = 0, of period 2ir 
(as expected); for A(j> E (0, 7r), the wave is shifted to the right, while for A<p E (ir,2ir), it lies to the 
left of the origin. In this case, we examine both the prediction for the relevant unstable eigenvalue, 
as well as the prediction of our theoretical results for the location of the center of the structure. The 
numerical results once again align extremely well with the theoretical ones, confirming the validity of 
our theoretical findings. 

Finally, we also considered a case in the neighborhood of the complete mutual cancellation of the two 
contributions of the effective potential. In particular, for the case of k\ = k<i = 2tt/5, and for B = 1, 
A(j) = 0, we examined the dynamics for A ~ 1.33, whereby in accordance with Fig. [H the relevant 
(translational) eigenvalue is numerically found to be returning to A 2 = 0, thereby restoring a regime 
of effective translational invariance. In this setting, we initialized a stationary soliton boosted by a 
factor of exp(iKx), with k = 0.1, in Fig. 0] (other values of k were also used with similar results). We 
observe that the soliton appears to propagate with a speed near the originally "assigned" one, being 
submitted only to very weak modulations due to the very weak (in this case) effective potential. These 
modulations are accompanied by oscillations of the solitary wave amplitude and lead to a speed slightly 
larger than 0.1 (shown by thick dashed line in Fig. H]). We note that in this way we can induce the 
robust motion of the waves over the combined linear and nonlinear lattice terrain. 

4 Conclusions 

In this paper, we have examined the evolution of bright solitary waves in the presence (and competing 
effects) of linear and nonlinear lattices. We have computed the effective potential landscape that the 
wave encounters and have explained how its curvature is associated with the wave stability. The 
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Figure 4: Space-time contour plot of the soliton evolution over a linear and nonlinear lattice with 
A = 1.33, B = 1, ki = &2 = 27r/5 and A(f> = 0. The thick dashed red line shows the curve x = O.lt for 
comparison. Notice the very weak modulation in the near-free propagation of the solitary wave. 



relevant translational eigenvalue has been explicitly computed and the transitions from stability to 
instability due to the zero crossing of this eigenvalue have been quantified as a function of the system's 
parameters. The same threshold condition has been obtained independently from a direct eigenvalue 
calculation. It has been shown that these theoretical frameworks capture accurately the location of the 
stationary waves, as well as the pertinent eigenvalues, hence they constitute valuable tools for inferring 
existence and stability information about the coherent structures of such models. On the other hand, 
the ensuing potential energy landscape can be used to derive a dynamical equation for the motion of 
the soliton; however, there is a number of notes of caution that should be made in that regard, as the 
latter may not capture entirely accurately the dynamical behavior, especially in "marginal" cases, due 
to the role of radiative corrections. It is also possible to appropriately tune the system parameters so 
as to nearly mutually cancel the effects of the linear and nonlinear lattice and produce a wave that is 
propagating at nearly-constant speed. 

It would be especially interesting to apply similar considerations to the case of dark solitons in 
repulsive BECs and examine their impact on the spectrum in the spirit of the recent work of [36J . Similar 
considerations could subsequently be extended to higher dimensional settings, where an important 
relevant example would be the influence of the lattice potentials to the existence and stability of 
structures with vorticity [371 EE] • 
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